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We investigate the entanglement properties of the joint state of a distinguished quantum system and its 
environment in the quantum Brownian motion model. This model is a frequent starting point for investigations 
of environment-induced superselection. Using recent methods from quantum information theory, we show that 
there exists a large class of initial states for which no entanglement will be created at all times between the 
system of salient interest and the environment. If the distinguished system has been initially prepared in a pure 
Gaussian state, then entanglement is created immediately, regardless of the temperature of the environment and 
the non-vanishing coupling. 
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No quantum system is completely isolated from its envi- 
ronment. This basic yet fundamental observation has been 
one of the key insights allowing an appropriate understanding 
of the dynamical emergence of classical properties in quan- 
tum systems. Not all initial states are equally fragile under 
the interaction of a distinguished quantum system with its en- 
vironment, and a relatively robust set of so-called preferred 
or pointer states is selected dynamically, a process typically 
referred to as environment-induced superselection (einselec- 
tion) or simply decoherence. This process is thought to play 
an important role in the transition from quantum to classical 
The most frequently employed model in investigations 
of einselection is the quantum Brownian motion model [HQ]. 
In this model one considers a distinguished quantum oscillator 
which is linearly coupled via the position operators to an en- 
vironment consisting of many harmonic oscillators. Initially, 
the state of the system of interest and its environment are as- 
sumed to be uncoiTelated, and the state of the environment is 
taken to be the canonical (Gibbs) state with respect to some 
temperature. The typical argument is that starting from the 
initial situation, the product state of the composite quantum 
system turns into a coiTelated state due to the interaction. If 
one considers the reduced state of the distinguished system 
one finds that it undergoes dissipation and decoherence. In 
the context of quantum Brownian motion it is often argued 
that entanglement is unavoidable. 

It is the aim of this letter to revisit the question of the 
creation of entanglement in quantum Brownian motion with 
recent powerful methods from quantum information theory 
[^,|[]. Our analysis will be split into two parts. In the first 
part we will show that surprisingly, quantum Brownian mo- 
tion does not necessarily create entanglement between the dis- 
tinguished system and its environment. The joint state of the 
system and its environment may be separable at all times, that 
is, not entangled [[7pO|]: All correlations are merely classical 
in the sense that one could prepare the same state by mixing 
product states, which can in turn be prepared by implementing 
local quantum operations only. By definition, separable states 
do not violate any Bell inequality. We then explicitly construct 
initial states with the property that no entanglement is created: 



they are mixed Gaussian states which are nevertheless differ- 
ent from Gibbs states. In contrast to the finite-dimensional set- 
ting, where a high degree of mixing automatically implies sep- 
arability [pi]], the existence of such initial states is not obvious 
l|T^]. The second part of our analysis is concerned with the 
question whether there exist initial states of the distinguished 
oscillator for which the joint state becomes immediately en- 
tangled. This question will be answered positively, and it will 
be demonstrated that all pure Gaussian states have this prop- 
erty, regardless of the initial temperature of the environment. 

From now on the distinguished quantum oscillator will be 
called S, the environment will be referred to as E. In the 
quantum Brownian motion model [^^, the total Hamiltonian 
consists of three parts H = H''^ ig) 1 + 1 ig) + H' , where 

2 ^+1 

2mi 2 ^ — ' 

h-=y:(-^p!+^xA. (1) 

The frequencies wi,... ,ujn+i and coupling constants 
K2, ■ • • , are taken to be positive. For convenience, we 
set uii = 1, all masses to be equal, rrij = 1 for j = 
1, . . . , + 1, and we require that the (A^ + 1) x (A^ + 1)- 
matrix V corresponding to the potential energy is positive, 
where Vi^i = 1^1/2 and Vj^j = uj'j/2, Vi^j = Vj^i = 
—Kj for j — 2, ...,N + 1, and all other entries of V are 
zero. Typically, one assumes product initial conditions [||Q], 
Po = Po ® Pa^ where the environment is initially in the 
Gibbs state p§ = exp(-/3_ff^)/tr[exp(-/3i?^)] associated 
with some inverse temperature (3. This model together with 
the above additional assumptions will be later referred to 
as QBM model in the more specific sense. The time evo- 
lution of the reduced state with respect to S can be deter- 
mined without approximations [Q]: for all spectral densities 
-^(^) = '^'j=^2^ KjS{uj — LUj)/{2ujj) one can derive a differ- 
ential equation that specifies the dynamical map. This com- 
pletely positive map £t, t S [0, oo), maps an initial state pif of 
S on the state pf = £t{Po) = tr£;[J7t(pf (g> Pq)U}] at a later 
time t, where Ut := sxp{—iHt). 
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We will first clarify the notation that will be used subse- 
quently. It will turn out to be appropriate not to investigate the 
state on the infinite-dimensional Hilbert space TL — TL^ ® TL^ 
of the joint system directly, but rather its associated covari- 
ance matrix. Throughout the paper we will make repeated 
use of the formalism of covariance matrices and their ma- 
nipulation by means of symplectic transformations [|l3 14 1. 
The 2?! canonical self-adjoint operators corresponding to po- 
sition and momentum of a system with n degrees of free- 
dom can be collected in a row vector O — {Oi , ■ ■ ■ , 02n) — 
{Xi ,Pi,... , Xn, Pn). The canonical commutation relations 
(CCR) can then be written in matrix form as [Oj,Ok\ = 
i(S2n)j,fc, giving rise to the skew-symmetric block diagonal 
real 2n x 2n-matrix Yi2n (or simply E when the size of the 
matrix is clear from the context). Gaussian states, which are 
defined through their property that the charateristic function 
is a Gaussian function in phase space, can be characterized 
in a convenient way through their moments. The first mo- 
ments {Oj)p, j — 1, . . . , 2n, are the expectation values of the 
canonical coordinates. The 2n x 2n covariance matrix F, 



r,- fc = 2tr [p (O, - (0,)p) {Ok - (Ofc)p)] - ^S„fe, (2) 

satisfying the Heisenberg uncertainty principle F + iS > 0, 
embodies the second moments of a state p. For later con- 
siderations we give the set of covariance matrices of a sys- 
tem with n degrees of freedom the name C2n := {F G 
M2n ■■ r = F^,r + zE > 0}, where M2„ is the set of 
real 2n x 27i-matrices. Particularly important will be co- 
variance matrices of Gibbs states cxp(— /3i7)/tr[exp(— 
with respect to a Hamiltonian H and the inverse tempera- 
ture f3, which are important examples of Gaussian states. For 
brevity, the corresponding covariance matrix will from now 
on be denoted as T{f3H). We will frequently employ linear 
transformations from one set of canonical coordinate to an- 
other which preserve the CCR, meaning that S'SS'^ — S. 
Such transformations form the group of (real linear) symplec- 
tic transformations 5*^(2^,11). A symplectic transformation 
S € Sp{2n, R) results in a change of the covariance matrix 
according to F i — > STS^, on the level of the states it is asso- 
ciated with a unitary operation p i — > U{S)pU (5*)^. 

We are now in the position to state the first Proposition. It 
is concerned with the fact that there exist initial states of the 
distinguished system and an initial temperature of the envi- 
ronment such that the joint state stays separable for all times. 
In the second part of the proof a lower bound for the required 
inverse temperature of the bath will be determined. 

Proposition 1. In the QBM model, for any choice of cou- 
pling constants (k2, . . • , kat+iJ, Kj > 0, and any frequencies 
{uj2, . . . , lun+i), > 0, there exists a Gaussian initial state 
of the system Pq with covariance matrix Fq and an inverse 
temperature f3 > of the environment such that 

Pt = Utipl ® exp{-pH^)/tr[exp{-pH'')])U} (3) 



is not entangled for all times t 6 [0, oo). Let ujoo ■ 
max{w2,-- - ,ojn+i}, S := 2^ - "^^k^, and il. := (w^ 



2^1/2)1/2^ f/zen the above inverse temperature is bounded 
from below by the smallest [3 that satisfies Fq ® T{f3H^) > 
T{-/H), where 7 min{2,log(l + 2/n)/n}. 

Proof of Proposition 1. The system S and environ- 
ment E together form a system with + 1 canonical de- 
grees of freedom, where now O = (Oi,... ,02W+2) = 
{Xi,Pi,X2, P2, ■ • ■ , Pn+i)- There exists a T e Sp{2N + 
2, R) such that the Hamiltonian H in the new canonical coor- 
dinates O"'^ = TO^ is the Hamiltonian 



1 q2 



N+1 p.2 



(4) 



of iV + 1 uncoupled oscillators, with real numbers ^"j /2, j = 
1, . . . , 7V+ 1, which are the eigenvalues of the positive matrix 
V . Due to the fact that the coupling is restricted to the coordi- 
nates associated with positions, T is both orthogonal and sym- 
plectic, i.e., T £ Sp{2N + 2, R) r\0{2N + 2). This symplec- 
tic transformation can be achieved by applying an appropriate 
orthogonal matrix AI G SO{N + 1) on the canonical vari- 
ables corresponding to position, (Oi, O3, . . . , 02n+i)'^ = 
M{Oi, O3, . . . , 02^+1)"^, such that H becomes diagonal in 
position, and the same matrix M on the momentum vari- 
ables, (02,04, . . . , d2N+2V = M(02, O4, . . . , 02JV+2)^. 

In these canonical coordinates the covariance matrix corre- 
sponding to the Gibbs state w.r.t. S becomes 

r{SH) = diag(F(<5#)i 1, . . . ,r{SH)2N+2,2N+2), 
nSH),^-i,2j-i - f{S^s)/^j, n5H)^,a, = fi5^,)^,^ (5) 

i = 1, . . . , iV + 1, where the function / : R+ — > R+ is 
defined as f{x) = 1 + 2/(exp(x) — 1). Let 7 > be defined 
as above. We now show that 



F(7i7) > 127V+2. 



(6) 



The largest eigenvalue of V is given by its operator norm 
Ip^] ||y||. By adding and subtracting the same term one 



1/2} 



V - 
I < 



obtains \\V\\ < max{Li;^/2, a;|/2, . . . , u./^_^ 
diag(w?/2,u;|/2, . . . ,cj^_^i/2)||, giving rise to \\V _ 
tj^/2 + (5^/^. The value of is related to the largest 
frequency cjoo := ma,x{LUj,j — 1, . . . , + 1} by \\V\\ — 
a)^/2. On using Eq. (Q), one finds after a few steps that in- 
deed f(7if) > I2W+2. 

In order to proceed, we need to invoke the concept of partial 
transposition. It has been shown in Ref. | ]T3| ] that in a system 
consisting of one oscillator in a system S and N oscillators in 
an environment E, a Gaussian state is separable if and only if 
its partial transpose is a quantum state. By using the matrix 
:— Tj2 ® (— S2Ar), the criterion can also be written in the 
form that a state with covariance matrix F is separable if and 
only if F + iS^^ > [O]. The next step is to see that with 
7 > as defined in Eq. (H 



T'^r(7i?)T + iE^^ > 127V+2 + iE'^^ > 0, 



(7) 
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since IliE^^II = 1 and T £ SO(2N + 2). Equipped 
with these preparatory tools, one can construct a class of 
product initial states such that the joint state of the system 
and its environment is separable at all times. These prod- 
uct states are of the form states p^, where = 
exp(— /3i?^)/tr[exp(— /3i7^)] with respect to a certain inverse 
temperature (3 > 0. p^ is already a Gaussian state, and the 
state of the distinguished system pg is taken to be Gaussian as 
well. On the level of covariance matrices such an initial state 
is represented as Tq = Tq © V{f-iH^), that is. 

To = © diag((ro)3,3, ■ • ■ , (ro)27V+2,2W+2), 
(ro)2i-l,2i-l = /(/3Wj)/Wj, (ro)2i,2i = f{l3ujj)ujj, (8) 

for j = 2,... ,N + 1, where £ C2. The covariance 
matrix Fq e C2 and /3 > are now chosen to be such 
that Fq - T'^T{'jH)T > holds. Such covai-iance matri- 
ces and inverse temperatures always exist. For a given co- 
variance matrix T^r{jH)T one can always choose (3 > 
and Ff e C2 in such a way that Fq - T'^T{'jH)T is diago- 
nally dominant JT^]. Since it is a symmetric matrix, diagonal 
dominance implies that Fq — T'^r{'yH)T is a positive ma- 
trix. From the definition of the covariance matrix r{jH) one 
can infer that Fo - T^r{-/H)T > is equivalent to the re- 
quirement F^ © T{[3H^) > F(7F), which is the inequality 
in Proposition 1 giving rise to the lower bound for /3. 

Since the Hamiltonian is a quadratic polynomial in the 
canonical coordinates, time evolution po 1 — > pt — UtPoUj is 
effected by a symplectic transformation on the level of covari- 
ance matrices. There exists a continuous map t e [0, 00) 1 — > 
St £ Sp{2N + 2, R), such that given a covariance matrix 
Fo £ C2JV+2 at time t — 0, the covariance matrix of pt be- 
comes Ft :— StTaSf . The aim is to show that this matrix 
corresponds to a separable state, i.e., Fj + i'S'^'^ > for all 
times. As r{'jH) is the covariance matrix of a Gibbs state 
with respect to H, StT{jH)S^ = T{-fH) for all t £ [0, 00), 
and hence, StT{-fH)Sj + iE^^ > 0. Moreover, the matrix 
5t(Fo - T{-fH))Sj' is positive, because Fq - T{-fH) IS posi- 
tive. It follows that 

Ft + iE^^ = StiTo - Ti-fH))Sj + T{jH) + lE^^ > 

for all t £ [0, 00). As the state was initially a Gaussian state, 
it remains Gaussian under time evolution. Having a positive 
partial transpose is equivalent with being separable for sys- 
tems where one of the parts consists of only one oscillator, 
which means that we can conclude that Ft corresponds to a 
separable state for all times. ■ 
So we have shown that for these initial conditions, no en- 
tanglement will be created at all times. At this point, a remark 
might be appropriate concerning the inverse temperature /3 in 
Proposition 1 . The question of the behaviour of (3 is partic- 
ularly relevant when one performs a continuum limit as is 
typically done when deriving quantum master equations. In 
this context, it is of interest to see that a lower bound for f3 
can be found that is independent on the number of oscillators. 
We consider a sequence of Hamiltonians {Hn}^^-^ of a joint 



system with an environment consisting of N oscillators, each 



equipped with a coupling constant and a frequency 
j = 2,...,iV + l. Take for each N an equidistant distribution 
of frequencies, such that ujj' — {j — l)ujao/N , where lOoo is 
the largest (cut-off) frequency. Concerning the spectral densi- 
ties we only make the assumptions that — aN^oj-^'Y with 
p > and a AT > 0, which covers the Ohmic {p = 1), the sub- 
ohmic {p < 1) and the supraohmic case {p > 1). The refine- 
ment must in all instances be made such that '}2!j=2 = 
2 J dujl{uj)uj remains constant. Then one can show that there 
exists a strictly positive lower bound for (3 as in Proposition 1 
which is independent of the number of oscillators N, and the 
limit iV ^ 00 may be performed. A detailed sketch of the 
argument will be given in footnote [[l^] . 

In turn, having this observation in mind one may ask 
whether there are initial states pg for which one can be sure 
that entanglement will be created immediately, no matter how 
weak the interaction is between the system and its environ- 
ment, and given any possibly very high initial temperature of 
the environment. We shall see that there exist such states: all 
pure Gaussian states have this property. 

Proposition 2. In the QBM model, for any initial 
pure Gaussian state pg of S, any coupling constants 
(k2, . . . , kn+i), Kj > 0, any frequencies {tui, . . . , cjat+i), 
ujj > 0, and any f3 > 0, the state pt — Ut{pQ (8" 
exp{—(3H^)/tr[exp{~(3H^)])Uf is entangled for all times 
t £ (0, e\for an appropriate e > 0. 

Proof. A Gaussian state pg with a covariance matrix 
Ff e C2 is pure iff (EsF)^ = -I2. The task is to find an 
e > such that £:t(p5'®exp(-/3i?^)/tr[exp(-/3iJ^)]) is en- 
tangled for all t £ (0,e]. Assume w.l.o.g. that {Oj)pg — 
for j — 1,2, i.e., all first moments vanish initially. Since 
{d/dt)\t=o{Oj) = 0, at t = the covariance matrix Ft of 
Pt = UtpoUf satisfies 



d_ 

dt 



t=o 



W^ToW 



(9) 



with some W £ M2N+2, which can be explicitly evaluated 
bymakinguseof itr[OjOfc(dpt/dt)|t=o] = ^[OjOk[H, po]], 
j, k = 1, . . . , 2N+2. The state pt is entangled at a time t > 
if the reduced state with respect to the system S and one of the 
oscillators of the environment E is an entangled state. The co- 
variance matrix of S and the oscillators of E with label 2 is 
given by a real 4x4 principal submatrix of Ft: its entries are 
(rt)j,fc with j,k — 1, . . . ,4. This matrix will from now on be 
called 7f e C4; it can be written in block form as 



7t 



At Ct 
CI Bt 



(10) 



where At,Bt,Ct £ M2- The statement that the reduced state 
of S and the oscillator with label 2 is entangled is equivalent to 
the statement that 7t +iE-^^ is not positive, where now E^« = 
E2 © (— E2). This in turn is equivalent with the smallest 
eigenvalue of (iE^^7t)^ being smaller than one (which means 



3 



that one of the symplectic eigenvalues of the partial transpose 
of 7t is smaller than one [p^). Let the smallest eigenvalue 
of (iEj^7t)^ be denoted as At, and let for brevity dt := 
det(At) +det(i?f ) — 2det(C(). The smallest eigenvalue Af can 
then be expressed as A* = (it/2-((i(/4-det(7f))^/^ jl6|. At 
t = 0, the covariance matrix 70 is of the form of Eq. (|10|) with 
Ao e C2, Co = 0, and Bq = diag(/(/3w2)/w2, /(/3t^2)w2), 
the latter matrix satisfying det(_Bo) > 1 by definition. As S is 
initially in a pure state, det(y4o) = 1. Therefore, Aq = 1. The 
first derivatives of At, Bt, and Ct at i = can be computed 
to be given by 



d_ 
di 
d_ 
di 
d_ 
dt 



At 



Ct = 



t=o 



Bt = 



(^0)2,2 - (^0)1,1 
(^0)2.2 - (Ao)i,i -2(Ao)i,2 

-K2(^o)l,l ' 
K2(^o)l.2 

(So)2,2 - (Boh.l 
(So)2,2 - (So)l,l 



(11) 



Hence, it follows that {dXt/dt)\t^Q < 0. This means that 
there exists an e > such that At < 1 for alH G (0, e], which 
implies that the partial transpose of the state corresponding to 
7t is not a state. We can conclude that the state pt associated 
with the covariance matrix Ff is an entangled state [ jl7| ] . ■ 

To summarize, we have investigated the entanglement prop- 
erties of the joint state of a distinguished system and its envi- 
ronment in quantum Brownian motion. Surprisingly indeed, 
we found that there exists a large set of initial states of the sys- 
tem for which no entanglement is created at all times. Also, 
we have shown that for pure initial Gaussian states of the dis- 
tinguished system entanglement will be immediately created. 
The tools we used were mostly taken from the field of quan- 
tum information theory. In fact, we hope that this letter can 
contribute to the line of thought of applying methods from 
quantum information theory to the issue of emerging classi- 
cality in quantum physics. 
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